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We adopt a field theoretical approach to the study of the structure and thermodynamics of a homogeneous 
Maier-Saupe nematogenic fluid interacting with anisotropic Yukawa potential. In the mean field approximation 
we retrieve a standard Maier-Saupe theory for liquid crystals. In this theory, the single-particle distribution 
function is expressed via the second order Legendre polynomial of molecule orientations. In the Gaussian 
approximation we obtain analytical expressions for correlation functions, free energy pressure, chemical po- 
tential, and elasticity constant. Subsequently we find corrections due to fluctuations and show that the single- 
particle distribution function now contains Legendre polynomials of higher orders. We also use Ward symmetry 
identities to set a simple condition for correlation functions. 
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1. Introduction 

Maier-Saupe nematogenic fluid [T] is one of the simplest models that account for the isotropic- 
nematic phase transition in the liquid crystal phase. The properties of this model have been inten- 
sively studied by the liquid theory methods such as integral equations for correlation functions . 
In the integral equation theory there is a problem of the correctness of taking the fluctuation effects 
into account, the treatment of which depends on closure relations used in integral equations. In 
order to treat the fluctuations more properly and to control the level of this treatment, in this paper 
we will apply the fleld theoretical approach. This is the flrst time the fleld theoretical approach is 
applied to the description of anisotropic molecular fluids. 

The method we are proposing focuses on fluctuations of the field at a given point and implements 
a perturbative scheme by expanding the Hamiltonian on density fluctuations. In the past, the 
statistical fleld theory proved to be successful in the description of a variety of systems with 
Coulomb and Yukawa-type interactions |13| fT^ . In this work we show that this approach 

also reproduces the familiar results for anisotropic systems, notably the mean fleld Maier-Saupe 
theory. Subsequently we go beyond this approximation and obtain an analytical expression for 
the pair correlation function. In the Gaussian approximation we also obtain new results for the 
main structural and thermodynamic properties of the system. The expressions we derive contain 
the orientational order parameter allowing us to compare the results for the isotropic and nematic 
phases. Finally, we calculate the correction to the mean fleld single-particle distribution function 
due to fluctuations which is expressed in terms of the fourth order Legendre polynomials of molecule 
orientations. Our results for the pair correlation functions predict the appearance of Goldstone 
modes in the system which is in full agreement with the theory of de Gennes [15j . 
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For the purpose of simplification, in this paper we consider a fluid of point particles. However, in 
the future we hope to modify the obtained results for non-point particles using the mean spherical 
results |5J|3J[3 as it was done for a non-point ionic system [16 . 

2. The model and field theory formalism 

We consider a molecular fluid of particles interacting via an anisotropic Yukawa-type potential 
iy{ri2 ,^1^12): 

y{n2.^i^2) = — e-"^i^P2(cos0i2) 

= — e-"'^-i^y2:„(i]i)r2„(i]2), (2.1) 

where ri2 denotes the distance between particles 1 and 2, 17 = (f?, (j)) are orientations of particles, 
P2(cos6'i2) = (3cos^ 012 — l)/2 is the second order Legendre polynomial of relative molecule ori- 
entations, Yira{^) are standard spherical harmonics ^17j without the normalization factor l/-v/47r, 
A is the amplitude of the interaction, and a is the inverse range. 

In a series of papers on ionic and Yukawa fluids [8] [13l [14] it was shown that it is possible to 
describe these fluids using the field theoretical approach. In this paper we will develop this approach 



for the description of an anisotropic molecular fluid with the interaction of the form (2.1 1. 

Within the field-theoretical formalism, the Hamiltonian is a functional of density field and can 
be written as 

/3i/[/7(r, n)] - /3i/™*'-[p(r, n)] + /?i7'"*[p(r, 17)] 

p(r, 17) [ln(p(r, f7)AR A|) - l] drdl7 

j v{ri2 , ^^if^2)p(ri , ni)p{r2 , f72)dridr2df7idl72 , (2.2) 



2 



where (3 = l/k^T is the inverse temperature, dfi = (l/47r) sin 6'd6'd0 is the normalized angle 
element, p{r,n) is particle density per angle such that J p(r, i7)df7 ~ p(r). At is the thermal de 
Broglie wavelength of the molecules, and the quantity is the rotational partition function for 
a single molecule fl7|. 

As in previous papers [SI [131 [14], we adopt the canonical ensemble approach. We fix the number 
of particles by the condition J p(r)dr = iV or ^ J p{r)dr = p, where V is the volume and p is 
the average density of the system. To verify this condition in a formally unconstrained calculus we 
introduce a Lagrange multiplier A such that 

The partition function Zj\j [/o(r, f7)] can be expressed as 

ZN[p{r,n)]= J Dp{r,n)exp{-l3H[p{r,n)]}, 

where Dp{v, il) denotes functional integration over all possible density distributions such that the 
total number of particles is N . The logarithm of the partition function gives the Helmholtz free 
energy 

l3F=-\nZN. (2.4) 

Due to the character of the interparticle interaction, the considered system is characterized by 
two non-dimensional parameters: non-dimensional density p* — pjo? and non-dimensional inverse 
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temperature /3* = —f3Aa = 1/T*. As we will see in our calculations, the third non-dimensional 
parameter appears M* = —AirpliA/a^ = Anp*f5* . 

In order to calculate the functional integral, we expand the Hamiltonian around the real angle- 
dependent density p{VL) which in the homogeneous case does not depend on r: 



l3H[p{Vt) + 5p{r, rj)] = 
= / {p{U) + 5p{v,n)) 



In 



In 1 



+ y{ri2.^i^2){p{^i) + 5p{v^,n^)){p{n2) + 5p{v2,n2))dv^dv2An^dn2. 



(2.5) 



3. Mean field approximation 



In order to obtain thermodynamic properties of the considered fluid we need to calculate the 
partition function. The lowest order approximation for the partition function is the saddle point 
for the functional integral which is the mean field approximation (MFA) from the physical point 
of view. In the canonical formalism it corresponds to fixing the Lagrange parameter A such that 
the relation (2.3 1 is true for the average density. 

we ol 

pm 



Expanding the logarithm in (2.5), we obtain 
5pH[p[v,^) 



= In- 



+ P j v{ri2 ,ni^2)p{^2)dv2dn2 . 



Sp{ri,n) 

The second term on the right-hand side of equation (3.1 1 equals 
13 J i^in2Mi^2)pi^2)dr2dfl2 = 



(3.1) 



y ^e^"'^-dr2 y ^^F2™'(^^2)1^2m'(f^i)E'''""^'™"(^2)da 

m' Im" 
Pl^^y^, P2mY2m{^l), 



where we have used 



A 

— ( 
ri2 



If we choose the value of parameter A to be 



AttA 



112 



J dQexp -(3^I^^p2rnY2rn{^) 



1 
Z 



then from (|3.l|) we get the following equation for density within MFA 

"/^F^y^P2m>2r«(^) 



P(^) = exp 



pfm, 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



where f{fl) is the single-particle distribution function and the averages can be calculated according 
to {...)n = J f{m---)dn. 

If we multiply both sides of equation (3.5 I by Y2mi^) and integrate by dfl we will obtain 

J^Pl'm' f dnYi,^,{n)Y2mm ^ p{Y2rrM)n . 
I'm' •' 

P2m ^ p{Y2m{^))Q ■ (3.6) 
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In normal ncniatics, the orientational distribution function /(fi) is axially symmetric with respect 
to a preferred direction n and depends only on the angle 6 between the molecular orientation n 
and n [2 . This means that only quantities independent of angle 0, in the plane perpendicular to n, 
yield non-zero averages and therefore for any m 7^ the averages {Y2rn{^))n equal 0. As a result, 
we obtain a well-known Maier-Saupe equation pL5] 



1 

Z 



Y2Q{rt) exp 



5 



an. 



(3.7) 



In terms of the orientational order parameter S and reduced unit M* , equations (3.5 1 and (3.7 1 
can be rewritten as follows: 



p{n) _ exp [^M*S cos^e] 

n " 1 

/ exp [|M*S'a;2] da; 



(3. 



S = {P2{cos0))n = ^{Y2o{n))n 



1 , 3i5 



/x2exp[|ArS'a;2] dx 



2 2 1 

/ exp [f Af*S'x2] da; 




(3.9) 



The order parameter S can take on values from to 1 with values S > corresponding to 



the nematic phase. Equation (3.9 1 is self-consistent and must be solved numerically. The resulting 
relationship between S and M* is presented in figure [T] The theory predicts a weak first-order 
phase transition from the isotropic phase with S" = to the nematic phase with S > 0. The 
smallest value of the order parameter S = 0.324 corresponding to M* — 4.484 defines the stability 
of the isotropic phase. A stable nematic phase is given by the solution that minimizes the free 
energy which in the MFA can be presented in the form 



N 



1 



= ln (pA^T^An) - I + -M* S 



(3.10) 



As a result, the stable nematic phase appears at M* = 4.587 and the value of the order parameter 
at the transition is 5* = 0.435. The region between M* = 4.484 and M* = 4.587 corresponds to the 
two-phase region which separates the isotropic and the nematic phases. The corresponding phase 
diagram in "density-temperature" coordinates is presented in figure [2] 
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4. Fluctuation and correlation effects: Gaussian approximation 

In the MFA, fluctuations are neglected. In this section we take them into account. To this end 



we should expand the Hamiltonian. From (2.5), the quadratic term in the Hamiltonian equals 

+ v{r^2 , n^n2)5p[v^ , n,)Sp{r2 , l]2)dridr2dr!idr!2 . (4.1) 
Expanding on the Fourier components 

^p(r,^!) = ^5p(k,^])e"'^ (4.2) 

k 

we obtain the expression for the quadratic term in the k-space 



4.1. Correlation functions 

The expression for a pair correlation function h is 

h{ri2 , f^if^2)(p(ri , f^i)) {p{Y2 , ^2)) - {5p{vi , ^ll)5p{v2 , ^2)) 

-(<Sp(ri , nO) {5p{v2 , ^2)) - 6{ri - r2)Sn,n, {p{vi , f^i)). (4.4) 

The second term on the right-hand side of equation ( |4.4[ ) disappears like in the homogeneous case 
{Sp{r, n)) — 0. The first term equals 

(Mri,f^i)Mr2,f^2)) = ^ / D{Sp{r,n)) e"^^^ W'-'f^)l<5p(ri , f^i)Mr2 , ^^2) 

J 

= 5^e*("-i-^^)((5p(k,r!i)5p(-k,172)) 

k 

^ /I}(,5p(k,f])) e-/3^^2[p(k,o)] ■ ^^-^^ 



As in the basis (4.2 1 the Hamiltonian is of diagonal form, and the Gaussian integral (4.5) yields 

(5p(k,l]0M-k,f^2)) = ^ (^ + ^^5E^2*™("l)^2™(1^2)) . (4.6) 

The inverse of the matrix in brackets is 

(^ + fc?f^^^^^"("^)^^"("^)) =Mfc,^^l^^2)p(^^l)p(^^2) + fe,o.p(^^l). (4.7) 



Identity ( |4.7[ ) is in essence the Ornstein-Zernike equation in the random phase approximation 
(RPA) for point particles [HI [H] 

h{k, niVi2) = c(fc, viiVL2) + j df73C(fc, nin3)h{k, ^^2^^3)/o(^^3) (4.8) 



33605-5 



M. Holovko, D. di Caprio, I. Kravtsiv 



with the closure 



C(ri2,f7if^2) = -M'^12,f^lf^2), 



(4.9) 



where C(fc,r2ir22) and /i(fc,r2ir22) are the Fourier transforms of the direct and pair correlation 
functions, respectively. 



In (4.8 1 we can expand f{k,iliQ2) on spherical harmonics 



(4.10) 



Due to the closure (4.9 1 and symmetry properties of the nematic we can write /(ri2 , in the 

form 



/(ri2 , ^1^2) = J2 f22miri2)Y^^i^i)Y2min2), 



(4.11) 



where 



C22m(r) 



1 

5 r 



This reduces to the following equation for harmonics 

h22m{k) = C22m{k) + (Kjm (^)) ^ P'i22m (fc) ^22™ (^) , 

resulting in the following expression for the harmonics of a pair correlation function 



1 



47r/3A 



h22m{k)- 5fc2+,2 + (y|^^(f^))^^l4,^^^. 

which is a renormalized, "effective" Yukawa potential in the k-space. In the r-space 

1 PA 



(4.12) 



(4.13) 



(4.14) 



h22m{r) 



5 r 



exp 



rxla2 + {Y2Un))nl47rpPA 



(4.15) 



where (KL(17))a = (1/p) / dl7p(f7)|r22„(f7)|. 

Dependencies of quantities (i2m(^))si the product SM* are presented in figure [s] Due to 
the normalization condition of functions Y2m{^); for SM*=0, the averages (i^2m(^))f2 = 1- 

4.2. Correction to the single-particle distribution function 

Correction to the single-particle distribution function due to Gaussian fluctuations can be found 
according to 



z 



(1 + [h{r, f7if}2) - C{r, nifl2)]) 



|M*Scos^ ( 



[l + Y.Y2Un)Y2mmh22m{r)+(3,yir)]) 



(4.16) 



where the normalization constant Z' can be found from condition J p'~^^ {il,)dft — p. Since 



\im[h22m{r)+Mr)]^~^P* 
1 — >-0 5 



l-^{Y2U^))nM* -1 



(4.17) 
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Figure 3. Dependence of quantities (l2m(f^))n on parameter for m = 0, 1, 2 

then the corrected single-particle distribution function has the form 



P 



/df7e 



M*S cos2 ( 



(4.18) 



We can also approximate the corrected single-particle distribution function in an exponential form 
as 



1 



Z' 



■ exp 



\m*S cos' ^ - E lY^mim' \ 1 - U^iMhM 



(4.19) 



where Z" is the normalization constant such that / p^-^^ {il)dn 
Note that 



(2/ + 1)1 





2 






2 




r 

\ m 


— TO 














F/o(f^), 



(4.20) 



where Z = 0, 2, 4; ( ^ q ) and ( o o o ) ^^'^ ^^-^^ corresponding Clebsch-Gordon coefficients fl7|. 

We can see that in the Gaussian approximation the dependence of the single-particle distribu- 
tion function on (3* and p* is more complicated than in the MFA: p{il) now depends not only on 
M* but there is also a direct /3*-dependence and a {Y2jji{^))n-dependence. We also see that in the 
Gaussian approximation the single-particle distribution function contains Legendre polynomials of 
the second and fourth orders of molecule orientations whereas in the linear approximation only 



Legendre polynomials of the second order are present. From expression (4.191 it is readily seen 



that the role of the fluctuation term increases with an increase of inverse temperature P* . 

4.3. Free energy, pressure, and chemical potential 

For a homogeneous system, the part of the Helmholtz free energy responsible for field interaction 
can be calculated by integrating with respect to the coupling parameter A: 

1 

F - i^id = ^ J drdrJidn2p(^^i)p(f^2)i'(r, flin2) J d\[l + /i(A, r, f^i^iz)], (4.21) 
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where 

ra 



E 



r 5 



exp 



r\a^ + {Yi^{n))n^^T:ppA\ 



Y;,^{n^)Y2^{n2). (4.22) 



Expression (4.211 in terms of parameter B — l^iirpPA yields 



HF~F,^)^V^{Y2om)l 



{a^ + {Y2Un))nBf' a{Y,Un))nB 
12tt Utt Sir 



Having an explicit expression for the free energy, we can find the pressure: 

T,N 



dV 



p+§^({Y2om)i+pi-^Y2om)i 



dp 



E 



{a^ + {Y,Un))nBf' a{Yi^{n))nB' 



where the derivatives of the averages (. . are equal to 



127r 127r 
d 



Stt 



d_ 
dp 



{Y2om)n 



B 


\Y20m) 


^,-{Y2l{n))^{Y2o{n)),, 


pa^ — pB 


\Y2om)i - {Y2%m)n 





In the isotropic phase {Y2o{ft))n = and expression (4.241 considerably simplifies: 



/3P = -/3 



T,N 



a^(a^ + B2m) ^'^ (a^ + B2m] 



1/2^ 



127r 



127r 



247r 



(4.23) 



(4.24) 



(4.25) 
(4.26) 



(4.27) 



where B2m = (^2m(^))o-S. This expression is similar to the one obtained in [131 [20], where there 
is a supplementary summation over m and the quantity B is replaced by i?2m • The chemical 
potential /i of the fluid can be found from expressions (4.23) and (4.24) as /i = (F + PV) /N and 
equals 

(3p - In (pAl Ar) + ^ (^{Y2om)l + \p§-p {Y2^m)j^ 

£ E (^(^2L(f^))a + ^^{YLm,^ (a - {a' + {YLmnPY 



/2 



(4.28) 



4.4. Broken symmetry problem and the elasticity constant 



A specific feature of the considered molecular fiuid is a broken symmetry which appears in the 
absence of an orienting external field. In [51|3J|7] using the Lovett-Mou-Buff-Wertheim equation |211 
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an exact relation for orientationally non-uniform fluids was obtained: 



Voi lnp{ni) - y"dri2dr!2C(ri2,r!ir!2)Vf 



(4.29) 



where C(ri2 ,il,iQ2) is the direct correlation function which in the RPA is given by equation ( |4.9[ ). 

Equation (4.29) is also known as the integro-differential form of the Ward identity [2, 7, 23]. 
The angular gradient operator Vq decomposes into 3 spherical components Vq , V_|_ , V_ . As Vq 
is oriented in the direction of liquid crystal, it vanishes due to rotational invariance. For other 
components, the following relations hold |17j : 



y±Yi^{n) = [l{l + 1) - m(m ± 1)]'/' Fi.^ii 
The direct correlation function can be expanded as j2] 



C(ri2,r!if^2) 



hnl'm' 



(4.30) 



(4.31) 



Due to the axial symmetry of a nematic, m = m' . From (4.29) we derive the following relation for 
the average {Y2i{^))n ■ 



C22ip(lfl(f^))i2 = 1, 



where 



C2 



dri2C22i(ri2) 



IAti(3A 



(4.32) 



(4.33) 



Due to condition (4.321, harmonic /i22i(fc) diverges in the limit k signalling the appearance of 



Goldstone modes in the system. This phenomenon is responsible for a number of unique properties 
of nematics such as elastic behavior and critical light scattering. The other harmonics of the pair 
correlation function should be finite according to the phenomenological theory of de Gennes [TS] . 
In the considered approximation (4.9) we have 



220 



C2 



C2 



and due to (4.32) and (4.33) 



1/2 



{Y2\m)n 



1/2 



(4.34) 



(4.35) 



As we can see from figure [3] for m = 2, the ratio (i^2^2(^))n/ {Yiii^))^ < 1 ^'^d h222{r) is well 
defined. However, for m = 0, the ratio (5^2o(^))o / (^2^i(^))n ^ ^ ^^'^ h22o{r) is not defined. For 
this reason, a problem arises also for thermodynamic properties and for the density correction. 
This problem is connected with approximation ( |4.9[ ) which leads to ( 4.34[ ). We should mention 
that in previous works [2J [31 [7] carried out in the mean spherical approximation equation (4.341 
is not true and the problem being discussed does not appear. We think that for point particles 
we can also solve this problem by introducing an additional isotropic interaction. We will consider 
this aspect of the problem in a separate paper. 

Broken symmetry is connected with elasticity properties of the considered fluid which are 
described by the elasticity constant K. Formal expression for K as proposed by Poniewiersky and 
Stecki )24| within our model reads 



/3K =^ I drdl7idr22?- 



2dpin,) dp{^l2] 



d cos 9i d cos 0; 



nx{^i)nx{^2)C{r, ^liU,2 



10 V 5' / r*C22i(r)dr 



(4.36) 
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Figure 4. Temperature dependence of the elas- 
ticity constant. 



Figure 5. Density dependence of the elasticity 
constant. 



Another way to calculate elastic constants comes from the theory of hydrodynamic fluctuations 
and can be found as 



f3K ^3 



lim fc^/i22i(fc) 



2 1 



(4.37) 



As we see, both expressions for K are identical. 

Dependencies of the reduced elasticity constant K* = —KjAc? = 127r5'^p*^ on temperature T* 
and density p* are presented in figures [4] and |5] One can see that with increasing temperature the 
elasticity constant decreases and with increasing density it increases. The effect is more pronounced 
respectively for larger density and for larger temperature. This result is in agreement with the 
previous result j7j obtained in the framework of the mean spherical approximation for a non-point 
model of the nematogenic Maier-Saupe fluid. 



5. Conclusions 

In this paper, for the flrst time the field theoretical approach is applied to the description of 
correlation functions and thermodynamic properties of molecular anisotropic fluids. As an example 
we consider a Maier-Saupe nematogenic fluid with the Yukawa potential of interparticle interaction. 
By expanding the Hamiltonian in powers of density fluctuations we examine the system in the mean 
field and Gaussian approximations. 

In the mean field approximation, we obtain analytical expressions for the single-particle distri- 
bution function and the orientational order parameter which are in agreement with the classical 
Maier-Saupe theory for nematic ordered fiuids. 

In the Gaussian approximation we find analytical expressions for the pair correlation function, 
the free energy, the pressure, the chemical potential, and the elasticity constant. We calculate the 
correction to the mean field single-particle distribution function due to fiuctuations and show that 
the corrected single-particle distribution function has a more complex dependence on density and 
temperature compared to the MFA. In contrast to the MFA single-particle distribution function, it 
includes the fourth order Legendre polynomials of molecule orientations in addition to the second 
order ones. We also use Ward symmetry identity to derive a simple expression for the average of 
spherical harmonic (Y2i{^))- One consequence of this condition is that for the system to be stable, 
the distance-dependent part of the interaction should be attractive. We show that in the Gaussian 
approximation the harmonic /i22i(^) diverges in the limit k — 0. Such a situation occurs at the 
phase transition from an isotropic to a nematic phase. This change in symmetry causes collective 
fluctuations known as the Goldstone modes. However, in the RPA for the considered system, the 
harmonic /i22o(^) is not defined. We hope to solve this problem in our next paper by introducing 
additional isotropic interparticle interactions. 
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HeMaTMHHMM n/iMH MaPiepa-Sayne: TeopeTMKO-no/ibOBi/iM niflxifl 

M. To/ioBKcP, P,. fli KanpicP, I. KpaB^i^ 

iHCTi/iTyr 4>i3MKii KOHfleHCOBaHi/ix CMCTeM HAH yKpaiHM, By;i. I. CBeHL4iL4bKoro, 1, 79011 JlbBis, yKpa'ma 

JlaBopaTopia e/ieKTpoxiMiT, xiMiT noBepxoHb i eHeprernHHoro MOfle/iioBaHHJi, 

BiflfliiieHHa xiMiT Bi/iLU,oT HaL4ioHa/ibHoV LUKO/11/1 riapiTex, n/i. >KyccK), 4, 75005 napi/i>K, 0paHL4i5i 

Mm sacTocoByeMO TeopeTHKO-nojibOBHW niflxifl p/\n Bi/iBHeHHa crpyKTypHi/ix i xepMOfli/iHaMinHi/ix BJiacTi/iBocxeM 
oflHopiflHoro HeMaxi/iHHoro n/inHy Mai/iepa-3ayne 3 aHisoTponHOK) BBaeMOflieio xi/iny KDKaBi/i. V Ha6/iM>KeHHi 
cepeflHboro nojia HaMi/i oxpwMaHO cxaHflapxHy xeopiio Mawepa-3ayne fl/i5i piflKi/ix Kpi/icxa/iiB. V xeopiT 
oflHOHacxwHKOBa cjiyHKLtin po3noflijiy BHpa>Kaexbca nepes noiiihOM /le>KaHflpa flpyroro nopiqflKy B3aeMHoT 
opieHxaL4iT nacxi/iHOK. V rayciscbKOMy Ha5jiM>KeHHi HaMM oxpuMaHO aHa/iixuHHi Bi/ipa3M p/\a KopejinL4ii/iHMX 
ctDyHKL4iM, KOHCxaHXM ejiacxMMHOcxi , Bi/ibHoT eHepriT, XHCKy i xiMiHHoro noxeHL4iajiy. 3a flonoMoroKi xoxo>KHOcxi 
Bopfla HaMM BcxaHOBJieHO npocxy yMOBy fljin Kope/iJiLiiwHi/ix cJjyHKuiw. HaMM xaKO>K 3HaHfleHO nonpasKH 
BHaciiiflOK 4>iiyKxyaL4iM i noKa3aHO mo BMpa3 fl/ia oflHOHacxMHKOBoV 4>yHKL4iT po3noflijiy B>Ke MicxMXb no^iiHOMH 
/le>KaHflpa bmlu,mx nopaflKiB. 

K/iiOHOBi c/iosa: HeMaTM4HMM nyiMH MaHepa-Sayne, reopeTiAKO-no/ibOBi/it^ niflxifl, Kopenau,iinHa (pyHKLtia, 
TepMOflUHaMiKa 
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